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1. Introduction 

n-point integrals in loop momentum space of tensor rank R appear in any realistic evaluation 
of Feynman diagrams. There are several ways to calculate them, and one approach expresses them 
by a small set of scalar integrals. The first systematic treatment, in a Standard Model calculation, 
is known as the Passarino-Veltman reduction ^ and expresses four-point tensor integrals (and 
simpler ones) algebraically by scalar one- to four-point functions. Two of us made use of this 
scheme in the early days [^, ^, |^. Nowadays tensor reductions became again a topic of research 
because at LEP2, LHC and ILC the interesting final states typically consist of more than two 
particles, some of them being massive. Several dedicated tensor reduction packages have been 
developed for the calculation of five- and six-point functions. As open-source packages we like to 
mention the Fortran packages LoopTools/FF ||5|, ^ (covering with n < 5,R < A) and Golem95 
(covering with n <6 and massless propagators) and the Mathematica package hexagon.m 
% |] (covering /,f with («,/?) < (6,4), (5,3)). 

Here, we describe a recursive implementation for tensor functions with n < 6,R < n for 
arbitrary internal masses. We use the Davydychev-Tarasov approach where the tensor integrals are 



first expressed by scalar integrals with higher dimensions and indices []10|]. In a second step, the 
scalar integrals may be expressed algebraically by scalar one- to four-point functions [ [TTl ] quite 
similar to the Passarino-Veltman reduction. In fact, when using the same basis, the approaches are 
equivalent. A difference, though, may arise in the algorithmic realization, and as a consequence in 
the numerical stability and speed of an implementation. For more comments on the differences of 



tensor reduction schemes we refer to the literature quoted. In a recent paper [12], we introduced a 
convenient and easy-to-program version of the reduction a la Davydychev-Tarasov, which allows a 
recursive determination of a chain of tensors. Here we will describe that scheme and present some 
numerical results. The integrals to be evaluated are: 

where the denominators cy have indices Vj and chords qy. 

Cj = {k-qjf -m] + ie. (1.2) 
The normalization C(£) plays a role for divergent integrals only and is conventional: 

Here, we use d = A — 2e and = 1 . For the evaluation of the scalar functions we will rely on either 
LoopTools or QCDloops/FF ||l3|, |^, and the latter one uses also C{e) as defined here. 

2. Recursions 



Our recursions begin with six-point functions where the well known formula [ ]14| , \\5\ |16[ 
may be used: 
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Figure 1: The recursion triangle. 



The auxiliary vectors Qs read: 



6 ('Oi\ 

e?=I^f7iP , ^=l--6. (2.2) 



'■=1 We 



The/j;_j' is obtained from /„ by shrinking line 5, and the (^y' J are signed minors of the 
modified Cay ley determinant ()^ [|T7|]. For further details of notations we refer to [12]. 

The further calculational chain may be read off from figure [I[ Its basic idea is to represent an 
«-point tensor by an «-point tensor of lower rank l^^^ and by all the {n — 1) -point tensors of 



lower rank ' For 5-point functions, we derived in [12|: 



^ /^'i-M^-.g^ - £/f •••^''^-''^eA', (2.3) 
= LliTT^ s = 0,...,n. (2.4) 

(=1 \)n 



The formula is the analogue to (2.1). For «-point functions with n <5, the corresponding represen- 
tations contain additional terms because the number of independent chords is then less than four 
so that the chords don't form a complete basis for d = 4. There are several modifications to be 
applied, and we hke to reproduce only one example with auxiliary terms: 

4 

t=\ 



' Similar recursive realizations of the Passarino-Veltman reduction may be found in |p^, see there figures 2 and 3. 
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P\ 


0.5 


0.0 


0.0 


0.5 


Pi 


0.5 


0.0 


0.0 


-0.5 


P3 


-0.19178191 


-0.12741180 


- 0.08262477 


-0.11713105 


PA 


-0.33662712 


0.06648281 


0.31893785 


0.08471424 


P5 


-0.21604814 


0.20363139 


-0.04415762 


-0.05710657 




P6 


= -(Pi +^2+^3 + ^4 + ^5) 





Table 1: Phase space point of massless six-point functions taken from 

with the additional tensor and vector components: 
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T^^ = (2o - L Qt - G^^ll > 
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^3 ~ ^3 tiO 2-1 2 ^« ' 

H=l 
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« = 0,...,4. 



i=\ 



(2.6) 
(2.7) 
(2.8) 
(2.9) 
(2.10) 



They may be, finally, represented by the scalar integrals in d dimensions 14,1^,12" ^^i ■> where the 
indices t,u,w indicate truncations of corresponding lines: 
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(2.11) 
(2.12) 
(2.13) 
(2.14) 



The representations for the simpler tensors have been given in [12]. 

3. Numerical results 

An example of a completely massless tensor reduction uses the momenta given in table |. 
We combined our tensor reduction with the scalar master integrals from QCDloop [13]. Table ^ 
contains sample tensor components of a six-point function with rank R = 5. It shows an agreement 
of eight digits between the results of our package Hexagon.F [ p^ and those of Golem95 for the 
constant terms of the tensor components. 
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Hexagon.F 


Golem95 


^03121 


0.158428987E+0, 0.41670698E-1 


0. 15842898 lE+0 , 0.41670700E-1 


^11020 


-0.143913860E+1, - 0.16464705E+0 


- 0. 1439 1 3 85 3E+ 1 , - 0. 1 6464708E+0 


^20200 


0.242928780E+2, 0.55504184E+2 


0.242928776E+2 , 0.555041 82E+2 


^22130 


0.225563941E+0, 0.23192857E+0 


0.225563949E+0 , 0.23192851E+0 


^33333 


0.244568 135E+0, 0.74014604E+0 


0.244568 138E+0 , 0.740146 lOE+0 



Table 2: Real and imaginary parts of selected tensor components of ranki? — 5 massless hexagon integrals; 
comparison of the packages Hexagon.F and Golem95. 



k + q. 



k + q2 



k + q4 




Figure 2: Momenta flow for the massive six-point topology. 



As a second example, we reproduce components of the massive tensor integral i^P^^^ in table 
Q The kinematics is defined by figure ^ (with = 0) and table ^. All tensor components are finite. 
The numbers could not be checked by another open source program, but we had the opportunity to 



compare them with an unpublished numerical package [19] 



Further numerical results may be found in the transparencies of the talk []20|]. 
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Pi 


0.21774554 E+03 


0.0 


0.0 


0.21774554 E+03 


P2 


0.21774554 E+03 


0.0 


0.0 


-0.21774554 E+03 


P3 


-0.20369415 E+03 


-0.47579512 E+02 


0.42126823 E+02 


0.84097181 E+02 


P4 


- 0.20907237 E+03 


0.55215961 E+02 


- 0.46692034 E+02 


- 0.90010087 E+02 


P5 


- 0.68463308 E+01 


0.53063195 E+01 


0.29698267 E+01 


-0.31456871 E+01 


P6 


-0.15878244 E+02 


-0.12942769 E+02 


0.15953850 E+01 


0.90585932 E+01 




mi = 


= m2 = m3 = mj = me 


= 110.0, m4= 140.0 





Table 3: Randomly chosen phase space point of six-point functions with massive particles. 





Hexagon.F 


^03121 


0.29834730E-09, - 0. 68229 122E-10 


y7 11020 


0.42830755E-09, 0.425748 11 E-09 


^20200 


-0.71172947E-08, 0.10102923E-07 


^22130 


- 0.29200434E-09, 0.785538 1 lE-10 


^33333 


0.17451484E-07, - 0.309 1431 6E-07 



Table 4: Selected tensor components of rank R ~5 massive hexagon integrals; 
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